First principles study of density, viscosity, and diffusion coefficients of liquid MgSiC>3 

at conditions of the Earth's deep mantle 
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Constant-pressure constant-temperature ab initio molecular dynamics simulations at high tem- 
peratures have been used to study MgSiC>3 liquid, the major constituent of the Earth's lower mantle 
to conditions of the Earth's core-mantle boundary (CMB). We have performed variable-cell ab initio 
molecular dynamic simulations at relevant thermodynamic conditions across one of the measured 
melting curves. The calculated equilibrium volumes and densities are compared with the simulations 
using an orthorhombic perovskite configuration under the same conditions. For molten MgSiOs, 
we have determined the diffusion coefficients and shear viscosities at different thermodynamic con- 
ditions. Our results provide new constraints on the properties of molten MgSiOs at conditions 
near the core-mantle boundary. The volume change on fusion is positive throughout the pressure- 
temperature conditions examined and ranges from 5% at 88 GPa and 3500 K to 2.9% at 120 GPa 
and 5000 K. Nevertheless, neutral or negatively buoyant melts from (Mg,Fe)Si03 perovskite com- 
positions at deep lower mantle conditions are consistent with existing experimental constraints on 
solid-liquid partition coefficients for Fe. Our simulations indicate that MgSiOs is liquid at 120 GPa 
and 4500 K, consistent with the lower range of experimental melting curves for this material. Linear 
extrapolation of our results indicates that the densities of liquid and solid perovskite MgSiOs will 
become equal near 180 GPa. 
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1. INTRODUCTION 

Melting is a ubiquitous process in planetary interiors 
and one of the dominant mechanisms for thermal trans- 
port and chemical differentiation in planets. The prop- 
erties of silicate liquids are thus essential for understand- 
ing a wide range of geophysical phenomena related to the 
deep Earth and its origin and evolution. (Mg, Fe)Si03 
perovskite is the most abundant mineral in the Earth's 
deep mantle. Here we report the first ab initio simu- 
lation of the properties of liquid MgSiOs at conditions 
corresponding to the deep interior of the Earth. 

There are a number of lines of evidence that strongly 
indicate the Earth (and other terrestrial planets) were 
partially or wholly molten at least at certain intervals 
during the accretion process 0,0- The subsequent cool- 
ing and crystallization of a magma ocean may lead to 
chemical differentiation of the mantle 3]. Thus, under- 
standing the dynamics of a terrestrial magma ocean is 
essential to understanding the initial conditions for the 
thermal and chemical evolution of the Earth. The liq- 
uid viscosity and its depth dependence are among the 
important parameters that characterize an early (liquid- 
rich) magma ocean and melt/crystal density inver- 
sions could greatly modify the structure and cooling his- 
tory of such an ocean 0. At present, only low-pressure 
experimental data or calculations based on empirical po- 
tentials jj] have been used to constrain these models. 

The existence of small degrees (< 1%) of partial melt 



in the present lower mantle has been suggested as an 
explanation for the strong relative variations of shear ve- 
locity relative to compressional (dlnVs /dlnVp) velocity 
in the deep mantle |5j. In the D" region at the base of 
the mantle, a seismic ultra-low velocity zone (ULVZ) of 
thickness 5-40 km has been detected locally on top of 
the core- mantle boundary 0, These regions are 
characterized by seismic compressional and shear veloc- 
ity reductions of ~ 10% and ~ 30% respectively. The 
presence of relatively large degrees of partial melt (~ 5- 
30%) has been proposed as the most plausible explana- 
tion for these features 0. There is also evidence locally 
for smaller shear velocity reductions in D" plausibly con- 
sistent with lesser amounts of partial melt at depths as 
great as 300 km above the core mantle boundary 0] . 

The presence of melt in the deep mantle can greatly 
affect a number of physical properties of the region. The 
viscosity of the melt may strongly modify heat transport 
and convective circulation within the boundary layer [jj . 
The presence of a partially molten layer suggests the den- 
sity contrast between solid and melt should be small. The 
buoyancy of the melt will be controlled by the intrinsic 
density difference, as well as compositional differences 
(i.e., Fe enrichment) between the melt and solid. Density 
inversions between silicate melts and equilibrium liquidus 
crystals have been extensively studied at upper mantle 
conditions 0,^3- However, there are only limited con- 
straints on possible density inversions under conditions 
of the Earth's lower mantle. Recently, it has been pro- 
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posed based on laboratory measurements to 15 GPa that 
basaltic melts may become denser than mantle peridotite 
at conditions near the base of the mantle However, 
these results are subject to considerable uncertainty due 
to the long extrapolations involved. More direct determi- 
nation of physical properties of the components of deep 
mantle melts are needed. 

Recent developments in the atomistic simulation of 
solids and liquids based on the full solution of the quan- 
tum mechanical equations for the electrons allow the 
theoretical study from first principles, i.e. without em- 
pirical or adjustable parameters, of the structural, ther- 
mal, and elastic properties of minerals at arbitrary con- 
ditions of pressure and temperature. Recent works have 
addressed successfully the thermoelastic properties of 
solid M gSiOa p erovskite at Earth's mantle conditions 
[Til [T5I Il6t flTl ITsI ]. Here we focus on the properties 
of liquid MgSi03 at high pressures, which we study us- 
ing first-principles molecular dynamics at constant pres- 
sure. We determine liquid densities and compare them 
with solid perovskite densities at similar conditions to 
extract melting properties. We also calculate or estimate 
dynamical properties of liquid MgSiC>3 such as diffusion 
and viscosity. 



2. TECHNICAL DETAILS 

Ab initio molecular dynamic simulations |l9ll are per- 
formed using density functional theo ry poll2l| within the 
generalized gradient approximation [22( (GGA). We use 
ultrasoft pseudopotential [2^ for O and norm-conserving 
pseudopotentials [24| for Mg and Si. Nonlinear core cor- 
rections [2^] were used for Mg. Kohn-Sham orbitals were 
expanded in plane waves with a kinetic energy cut off 
of 30 Ry. The simulation cell contained 80 atoms (16 
MgSi03 units) and a time step of 13 a.t.u. was used. 
Variable-cell dynamics were used to impose the required 
value of the pressure l2a[ . Temperature was controlled 
by a Nose thermostat |27l |28| . The liquid configuration 
was obtained by heating the system to 10000 K and then 
decreasing the temperature to 5000 K in a total time of 
1 ps. 

The diffusivity, or diffusion coefficient (D) is estimated 
by calculating the mean square displacement of the atoms 
as a function of time. The diffusion coefficient is given 
by the Einstein relation 29] : 



D« lim I(|r(t +t)-r(i )| 2 ) 

t->oo oTj 



lim ^B?{t)). (1) 



The calculated mean square displacement were fitted by 
a straight line: (5R 2 (t)) « 6Dt + b. The viscosity (77) 
is estimated using the generalized Debye-Stokes-Einstein 
formula derived by Zwanzig |30j | : 



TABLE I: Calculated properties of MgSi0 3 at 3500 K as 
a function of pressure using solid perovskite and liquid con- 
figurations. For liquid configuration at 88 GPa, atoms are 
diffusive and the estimated diffusion coefficient is 1.14 x 
10~ 5 cm 2 /s. 



Solid perovskite 
T P V s P s 
(K) (GPa) (A 3 ) (g/cm 3 ) 



Vl 

(A 3 ) 



Liquid 

PL AV 

(g/cm 3 ) (%) 



3500 88 

3500 93 

3500 100 

3500 110 

3500 120 

3500 135 



140.787 
139.076 
137.263 
134.786 
132.013 
128.815 



4.718 
4.776 
4.839 
4.928 
5.031 
5.156 



147.764 
145.296 
142.098 
138.681 
135.360 
131.762 



4.495 
4.571 
4.674 
4.790 
4.907 
5.041 



4.96 
4.47 
3.52 
2.89 
2.54 
2.29 



Dr] 



ksTn 1 / 3 



= 0.0658(2 + 77/771) = C/. 



(2) 



Here n is the number density of the atoms, n = N/V. 
The constant Ci depends on the ratio of the shear vis- 
cosity (rf) to the longitudinal viscosity (r/i). Although 
these numbers are not actually available, Ci has bounds 
that can vary between 0.132 and 0.181. In our work, we 
adopt a typical value of Ci = 0.171. Wasserman et al Q 
studied the transport properties of perovskite melts by 
molecular dynamics under temperatures (3500-6000 K) 
and pressures (5.1-78 GPa) and their results agreed rea- 
sonably well with Zwanzig's formula. The general va- 
lidity of Zwanzig's formula has been studied in detail in 
other systems |3l|, |22, |33 ■ 



3. SIMULATION RESULTS 

3.1. Materials properties at high pressures and 
temperatures 

The calculated equilibrium volumes and densities of 
four MgSiC>3 units at 3500 K as a function of pressure are 
given m Tabled The results are shown in the right panels 
of Fig. n Our results at 88 GPa with the liquid configura- 
tion showed that atoms are diffusive with a diffusion coef- 
ficient of 1.14 x 10~ 5 cm 2 /s. Even though this thermody- 
namic point (3500 K and 88 GPa) is close to the melting 
line obtained by the Berkeley group 0, |3^, |3(j , the sta- 
bility of the liquid in our simulations does not allow us 
to conclude that 3500 K is above the melting line of the 
theoretical model employed. Liquid (mcta)stability in 
the few picoseconds of a first-principles simulation could 
be due to super-cooling effects |2^| . We did not estimate 
the diffusion coefficients at higher pressures because all 
the available melting curves suggest that MgSi03 can 
only be solid at these thermodynamic conditions. 

Regarding the simulations with liquid configuration, 
when pressure increases, the atoms approach each other 
and become no longer diffusive. This can be observed by 
our estimated diffusion coefficients at 88 GPa (1.14x 10 -5 
cm 2 /s) and 120 GPa (0.04 x 10~ 5 cm 2 /s). In other words, 
the under-cooled liquid is glassified. This agrees with the 
melting curve obtained from Knittle and Jeanloz |35| . 
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FIG. 1: (Left panels) Volumes and densities of solid and 
molten MgSiOg at different temperatures, pressure is fixed 
at 120 GPa. Note the difference between T > 4500 K and 
T < 4000 K. (Right panels) The same quantities but with 
temperature fixed at 3500 K and pressure changes. There 
is a tendency of density inversion between melting solid and 
liquid at higher pressures. 



and Heniz and Jeanloz |36j |. In addition, as pressure 
increases, the percentage difference of volumes between 
solid and liquid MgSiOa becomes smaller, which indicate 
a tendency of density inversion at higher pressures. In 
the left panels of Fig. ^ the volumes and densities are 
shown as a function of temperature up to 5000 K. Pres- 
sure is fixed at 120 GPa. This is close to that at the 
top of the D" region in the mantle. This pressure is also 
close to where MgSiOa transforms to a post-perovskite 
phase [3J| In this study, our solid state calculations are 
restricted to the perovskite crystal structure. The volume 
of the post-perovskite (CalrOs-type) phase is ~1-1.5% 
less than perovskite jH, |3^, 0] a t deep lower mantle 
pressures, so this transformation will increase the density 
contrast between solid and liquid. Because of uncertain- 
ties associated with the Clapeyron slope of the transition 
and the deep mantle geotherm, it is possible that at the 
base of the mantle the geotherm will cross back into the 
perovskite stability field |4l|. 

Although melting has not been observed during the 
simulations, our results provide an upper bound of the 
solid volume. The numerical results are given in Table 1TT1 
For liquid simulations, the changes of volume and density 
can be approximately divided into two regions, say, T > 
4500 K and T < 4000 K. It may be due to the possible 
glassification of liquid below 4000 K. Similar observation 
of glassification can be found if one looks at the mean 
square displacements of the liquid at different tempera- 
tures. 
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FIG. 2: Mean square displacements of molten MgSiOs at 120 
GPa. The liquid starts to glassify at temperature around 4000 
K. 
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FIG. 3: Mean square displacements of molten (upper panel) 
and solid (lower panel) MgSi0 3 at 5000 K and 120 GPa. Solid 
lines: O atoms, dashed lines: Mg atoms, long dashed lines: 
Si atoms, solid bold lines: all atoms. 



3.2. Diffusion near the core-mantle boundary 

In Table ^ we tabulate the equilibrium volumes and 
densities of solid and molten MgSiC>3, together with 
the estimated diffusion coefficients and viscosities of the 
molten MgSiC>3 as a function of temperature. Pressure is 
fixed at 120 GPa. The volumes and densities are shown 
in the left panels of Fig. 2] The estimated diffusion co- 
efficient drops drastically as the temperature drops from 
4500 K to 4000 K. This is a strong indication of glassifica- 
tion of molten MgSiO*3. The mean square displacements 
of molten MgSiC>3 at different temperatures are shown in 
Fig. [2 At higher temperatures (4500 K ~ 5000 K), The 
atoms are clearly diffusive. This provides strong evidence 
that molten MgSi0 3 exists at 4500 K and 120 GPa. 

The mean square displacements of different species of 
atoms are shown in Fig. [3J Both solid and liquid results 
are shown. In the solid, Si atoms are bounded by the 
octahedra formed by the O atoms, and each octahedron 
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TABLE II: Calculated materials properties of MgSiOa at 120 GPa as a function of temperature using solid perovskite and 
liquid configurations. 







Solid 






Liquid 




T 


P 


Vs ps 


V L 


PL 


D 


V 


(K) 


(GPa) 


(A 3 ) (g/cm 3 ) 


(A 3 ) 


(g/cm 3 


) (lO^ 5 cm 2 /s) 


(CP) 


3500 


120 


132.013 5.031 


135.360 


4.907 


0.04 




4000 


120 


133.340 4.981 


135.712 


4.895 


0.46 


107.74 


4500 


120 


134.476 4.939 


138.448 


4.798 


1.77 


31.46 


5000 


120 


135.792 4.891 


139.795 


4.751 


3.33 


18.51 



is then surrounded by 8 Mg atoms. Hence Mg atoms are 
expected to have a larger vibrational amplitude then that 
of the other species, and Si atoms are expected to have 
the smallest vibrational amplitude (Fig. |3J. In molten 
MgSiOs, the slopes of the mean square displacements of 
each species is close to that of each other. This suggests 
the self-diffusion rate of each species is almost the same. 
The simplest explanation is that the three types of atoms 
diffuse like a single molecule, which has to be examined 
carefully. However, as the diffusion mechanism is beyond 
the scope of this work, it will be the subject of further 
studies. 



4. DISCUSSION AND CONCLUSION 

At 3500 K, the density of crystalline MgSiOs in the 
perovskite structure is 5% greater than the liquid at 88 
GPa, and this density excess is 2.9% at 120 GPa and 5000 
K. The transformation to a post-perovskite phase would 
further enhance this density contrast. Thus, we find that 
liquid MgSiOa is always less dense than the correspond- 
ing solid under conditions encountered in Earth' s man- 
tle. This implies that stability of melts in the deep lower 
mantle requires in addition compositional differences be- 
tween melt and solid. It is well known on the basis of 
mineral-melt partitioning data for perovskite that Fe is 
preferentially partitioned into the melt. Hence, in addi- 
tion to the volume of fusion, Fe enrichment in the melt 
must be considered when evaluating the stability of melts 
in an idealized (Mg,Fe)Si03 lower mantle [42| . 

The effect of Fe on densities in crystalline and liquid 
(Mg,Fe)Si03 perovskites can be calculated using known 
thermoelastic parameters 0> EH and assuming the ef- 
fect of Fc on density in the liquid is similar to its ef- 
fect on the solid. For solid (Mg,Fe)Si03 compositions 
with Mg/(Mg+Fe) ratios of 0.90 to 0.95, the solid-liquid 
Fe partition coefficients required to be in equilibrium 
with neutrally buoyant melt at 120 GPa and 5000 K 
range from 0.31 to 0.48. These are within the wide 
range of reported experimental partition coefficients for 
perovskite compositions at lower mantle pressures (0.12- 
0.57) [HI!! EH. Thus, the combination of the small, 
positive volume change of fusion determined here with 
Fe partitioning constraints from experiment suggest that 
negative or neutrally buoyant melts are plausible in the 
D" and ULVZ layers at least with respect to perovskite 



structure solids. Partition coefficients within the exper- 
imentally measured range (0.2-0.34) also produce neu- 
trally buoyant melts at 88 GPa and 3500 K. Thus, these 
results suggest that small amounts of neutrally buoy- 
ant partial melt in the deep lower mantle could also 
contribute to the anomalous <9 In Vs /<9 In Vp observed in 
this region |j. Negative buoyancy under these condi- 
tions could instead lead to melt pooling in D" and the 
ULVZ mi • Better constraints on Fe partition coefficients 
under deep mantle conditions are needed to distinguish 
between these possibilities. 

Extrapolation of the trend in Fig. 2] indicates that 
MgSi03 liquid and crystal densities will become equal 
near 180 GPa. There are few experimental constraints 
on densities in MgSiOs liquids under these extreme con- 
ditions. The density change along the Hugoniot in 
Mg2Si04 upon shock melting at 150 GPa is small and 
there is a suggestion of enhanced compressibility in the 
melt More recent shock experiments on MgSiOs 

compositions have suggested that the density of the melt 
becomes comparable to that of the solid near 120 GPa, 
and at 170 GPa the melt density exceeds the solid den- 
sity mi- However, the shock data near 120 GPa (which 
are directly comparable to our simulations) have uncer- 
tainties that allow melt densities to be as much as several 
percent less than the solid, and thus they are not neces- 
sarily inconsistent with the present results. 

The viscosities of silicate liquids have been investigated 
experimentally using a wide variety of techniques at am- 
bient and low pressures |4^. However, there are few 
constraints on the behavior of silica poor and relatively 
depolymerized melts that are relevant to melting in the 
deep Earth. The viscosity of a diopside (CaMgSi20g) liq- 
uid has recently been determined to high pressures and 
temperatures using the falling sphere method |50l | . At 8- 
13 GPa and 2200-2470 K, the reported viscosities range 
from 28-510 centipoise which are generally larger than 
our higher pressure and temperature values (Table [nj . 
The temperature dependence of viscosity for MgSiOs liq- 
uid remains strong even at 120 GPa (Table ITT|l . 

In this work, we preformed ab initio molecular dy- 
namic simulations to study MgSiOs, the main mineral 
composition in the Earth's lower mantle. By applying 
variable-cell dynamics, we studied both solid perovskite 
and molten MgSiOs in constant NPT ensemble without 
the need for optimizing the cell parameters. At ther- 
modynamic conditions close to that of the core-mantle 
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boundary, we have constrained the viscosity, diffusion co- 
efficient, and density of liquid MgSi0i3. Our results also 
support the validity of the relatively low melting curve 
for MgSiC>3 perovskite |34l l35l l36l| as compared with that 
reported by Zerr and Boehler [5l|. 

The next step is to study the diffusion mechanism in 
molten MgSiC>3 as well as the chemical heterogeneity on 
the diffusion and melting behavior. With the increasing 
computation power, it is anticipated that the applica- 
tions of ab initio molecular dynamics can contribute to 



detailed understanding of mineral physics at the core- 
mantle boundary. 



Acknowledgments 

J.T.K. Wan acknowledges the financial support from 
The Croucher Foundation. This work was supported by 
the NSF. We thank A. Kubo for valuable discussion. 



[1] D. J. Stevenson, Formation and early evolution of the 
Earth, in Mantle Convection: Plate Tectonics and Global 
Dynamics, edited by W. R. Peltier, pp. 817-873, Gordon 
and Breach, New York, 1989. 

[2] V. S. Solomatov, Fluid dynamics of a terrestrial magma 
ocean, in Origin of the Earth and Moon, edited by R. M. 
Canup and K. Righter, pp. 323-338, U. of Arizona Press, 
Tucson, 2000. 

[3] E. Ohtani, The primordial terrestrial magma ocean and 
its implication for stratification of the mantle, Earth 
Planet. Sci. Lett. 38, 70-80, 1985. 

[4] E. A. Wasserman, D. A. Yuen, and J. R. Rustad, Molec- 
ular dynamics study of the transport properties of per- 
ovskite melts under high temperature and pressure con- 
ditions, Earth Planet. Sci. Lett. 114, 373-384, 1993. 

[5] T. S. Duffy and T. J. Ahrens, Lateral variation in lower 
mantle seismic velocities, in High-Pressure Research in 
Mineral Physics: Applications to Earth and Planetary 
Sciences, edited by Y. Syono and M. H. Manghnani, pp. 
197-206, Terra Scientific, Tokyo, 1992. 

[6] E. J. Garenro, S. P. Grand, and D. V. Helmberger, Low 
p-wave velocity at the base of the mantle, Geophys. Res. 
Lett. 20, 1843-1846, 1993. 

[7] J. Revenaugh and R. Meyer, Seismic evidence of partial 
melt within a possibly ubiquitous low-velocity layer at 
the base of the mantle, Science 277, 670-673, 1997. 

[8] D. Helmberger, S. D. Ni, and J. L. X. Wen Ritsema, Seis- 
mic evidence for ultralow-velocity zones beneath africa 
and eastern atlantic, J. Geophys. Res. 105, 23865-23878, 
2000. 

[9] Q. Williams and E. J. Garnero, Seismic evidence for 
partial melt at the base of Earth's mantle, Science 273, 
1528-1530, 1996. 
[10] L. X. Wen, P. Silver, D. James, and R. Kuehnel, Seismic 
evidence for a thermo-chemical boundary at the base of 
the Earth's mantle, Earth Planet. Sci. Lett. 189, 141- 
153, 2001. 

[11] C. B. Agee, Crystal-liquid density inversions in terrestrial 
and lunar magmas, Phys. Earth. Planet. Int. 107, 63-74, 
1998. 

[12] S. M. Rigden, T. J. Ahrens, and E. M. Stolper, Densities 
of liquid silicates at high pressures, Science 226, 1071- 
1074, 1984. 

[13] E. Ohtani and M. Maeda, Density of basaltic melt at 
high pressure and stability of melt at the base of the 
lower mantle, Earth Planet. Sci. Lett. 193, 69-75, 2001. 

[14] B. B. Karki et al., Elastic properties of orthorhombic 
MgSiOa perovskite at lower mantle pressures, Am. Min- 



eral. 82, 635-638, 1997. 
[15] B. B. Karki, R. M. Wentzcovitch, S. de Vironcoli, and 
S. Baroni, Ab initio lattice dynamics of MgSiOa per- 
ovskite at high pressure, Phys. Rev. B 62, 14750-14756, 
2000. 

[16] A. R. Oganov, J. P. Brodholt, and G. D. Price, Ab 
initio elasticity and thermal equation of state of MgSiOa 
perovskite, Earth Planet. Sci. Lett. 184, 555-560, 2001. 

[17] A. R. Oganov, J. P. Brodholt, and G. D. Price, The 
elastic constants of MgSiOa perovskite at pressures and 
temperatures of the Earth's mantle, Nature 411, 934- 
937, 2001. 

[18] R. M. Wentzcovitch, B. B. Karki, M. Cococcioni, and 
S. de Vironcoli, Thermoelasticity of MgSiOs perovskite: 
insights on the nature of the Earth's lower mantle, Phys. 
Rev. Lett. 92, 018501, 2004. 

[19] R. Car and M. Parrinello, Unified approach for molecular 
dynamics and density functional theory, Phys. Rev. Lett. 
55, 2471-2474, 1985. 

[20] P. Hohenberg and W. Kohn, Inhomogeneous electron 
gas, Phys. Rev. 136, B864-B871, 1964. 

[21] W. Kohn and L. J. Sham, Self-consistent equations in- 
cluding exchange and correlation effects, Phys. Rev. 140, 
A1133-A1138, 1965. 

[22] J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized 
gradient approximation made simple, Phys. Rev. Lett. 
77, 3865-3868, 1996. 

[23] D. Vanderbilt, Soft self-consistent pseudopotentials in a 
generalized eigenvalue formalism, Phys. Rev. B 41, 7892- 
7895, 1990. 

[24] N. Troullier and J. L. Martins, Efficient pseudopotentials 
for plane-wave calculations, Phys. Rev. B 43, 1993-2006, 
1991. 

[25] S. G. Louie, S. Froyen, and M. L. Cohen, Nonlinear ionic 
pseudopotentials in spin-density functional calculations, 
Phys. Rev. B 26, 1738-1742, 1982. 

[26] M. Parrinello and A. Rahman, Crystal structure and 
pair potentials: A molecular-dynamics study, Phys. Rev. 
Lett. 45, 1196-1199, 1980. 

[27] S. Nose, A molecular-dynamics method for simulations 
in the canonical ensemble, Mol. Phys. 52, 255-268, 1984. 

[28] W. G. Hoover, Canonical dynamics: Equilibrium phase- 
space distributions, Phys. Rev. A 31, 1695-1697, 1985. 

[29] M. P. Allen and D. J. Tildesley, Computer Simulation of 
Liquids, Oxford University Press, New York, 1987. 

[30] R. Zwanzig, On the relation between self-diffusion and 
viscosity of liquids, J. Chem. Phys. 79, 4507-4508, 1983. 

[31] N. H. March and M. P. Tosi, Generalized stokes-einstein 



6 



relation for liquid metals near freezing, Phys. Rev. E 60, 
2402-2403, 1999. 

[32] B. Bagchi, Relation between orientation correlation time 
and the self- diffusion coefficient of tagged probes in viscus 
liquids: A densitional theory analysis, J. Chem. Phys. 
115, 2207-2211, 2001. 

[33] J. D. Gezelter, E. Rabani, and B. J. Berne, Calculating 
the hopping rate for diffusion in molecular liquids: CS2, 
J. Chem. Phys. 110, 4333-3452, 1999. 

[34] R. Jeanloz and A. Kavner, Melting criteria and imag- 
ing spectroradiometry in laser-heated diamond-cell ex- 
periments, Phil. Trans. R. Soc. Lond. A 354, 1279-1305, 
1996. 

[35] E. Knittle and R. Jeanloz, Melting curve of (Mg,Fe)Si0 3 
perovskite to 96 GPa - evidence for a structural transition 
in lower mantle melts, Geophys. Res. Lett. 16, 421-424, 
1989. 

[36] D. L. Heinz and R. Jeanloz, Measurement of the melting 
curve of Mg .gFeo.iSi03 at lower mantle conditions and 
its geophysical implications, J. Geophys. Res. 92, 11437- 
11444, 1987. 

[37] M. Murakami, K. Hirose, K. Kawamura, N. Sata, and 

Y. Ohishi, Post-perovskite phase transition in MgSiOa, 

Science 304, 855-858, 2004. 
[38] A. Oganov and S. Ono, Theoretical and experimental 

evidence for a post-perovskite phase of MgSiOa in Earth's 

D" layer, Nature 430, 445-448, 2004. 
[39] T. Tsuchiya, J. Tsuchiya, K. Umemoto, and R. Wentz- 

covitch, Phase transition in MgSiOs perovskite in the 

Earth's lower mantle, Earth Planet. Sci. Lett. 224, 241- 

248, 2004. 

[40] S. Shieh, T. Duffy, A. Kubo, G. Shen, and 
V. Prakapenka, Synthesis of post-perovskite phase from 
a natural orthopyroxene, EOS Trans. AGU 85, 2004, Fall 
Meet. Suppl., Abstract MR23A-0182. 

[41] J. Hernlund, C. Thomas, and P. Tackley, Post-perovskite 
double-crossing, partial melting, and the thermal struc- 
ture of Earth's D" layer, EOS Trans. AGU 85, 2004, Fall 



Meet. Suppl., Abstract MR22A-02. 

[42] E. Knittle, The solid/liquid partitioning of major and 
radiogenic elements at lower-mantle pressures: Implica- 
tions for the core-mantle boundary region, in The Core 
Mantle Boundary Region, edited by M. Gurnis, M. E. 
Wysession, E. Knittle, and B. A. Buffett, pp. 119-130, 
American Geophysical Union, Washington D. C, 1998. 

[43] I. Jackson, Elasticity, composition, and temperature of 
the Earth's lower mantle: A reappraisal, Geophys. J. Int. 
134, 291-311, 1998. 

[44] F. Deschamps and J. Trampert, Mantle tomorgraphy 
and its relation to temperature and composition, Phys. 
Earth Planet. Int. 140, 277-291, 2003. 

[45] E. Ito and E. Takahashi, Melting of peridotite at up- 
permost lower mantle conditions, Nature 328, 514-517, 
1987. 

[46] E. A. McFarlane, M. J. Drake, and D. C. Rubie, El- 
ement partitioning between Mg-perovskite, Geochim. 
Cosmochim Acta 58, 5161-5172, 1994. 

[47] J. M. Brown, M. D. Furnish, and R. G. McQueen, Ther- 
modynamics for (Mg,Fe)2Si04 from the hugoniot, in 
High-Pressure Research in Mineral Physics, edited by 
M. H. Manghnani and Y. Syono, pp. 373-384, American 
Geophysical Union, Washington DC, 1987. 

[48] J. Akins, S. Luo, P. Asimow, and T. J. Ahrens, Shock- 
induced melting of MgSiOs perovskite and implications 
for melts in Earth's lowermost mantle, Geophys. Res. 
Lett. 31, 2004. 

[49] R. H. Doremus, Viscosity of silica, J. Appl. Phys. 92, 

7619-7629, 2002. 
[50] J. Reid et al., The viscosity of CaMgSi2 06 liquid at 

pressures up to 13 GPa, Phys. Earth Planet. Int. 139, 

45-54, 2003. 

[51] A. Zerr and R. Boehler, Melting of (Mg,Fe)Si0 3 - 
perovskite to 625 kilobars - indication of a high-melting 
temperature in the lower mantle, Science 262, 553-555, 
1993. 



